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The first measurement of the charged-particle multiplicity density at mid-rapidity in Pb-Pb col- 
lisions at a centre-of-mass energy per nucleon pair ^/snn = 2.76 TeV is presented. For an event 
sample corresponding to the most central 5% of the hadronic cross section the pseudo-rapidity den- 
sity of primary charged particles at mid-rapidity is 1584 ± 4 (stat.) ± 76 (sys.), which corresponds 
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sions at y^SNN = 0.2 TeV. This measurement provides the first experimental constraint for models 
of nucleus-nucleus collisions at LHC energies. 
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The theory of strong interactions, Quantum Chromo- 
Dynamics (QCD), predicts a phase transition at high 
temperature between hadronic matter, where quarks and 
gluons arc confined inside hadrons, and a deconfmcd 
state of matter, the Quark-Gluon Plasma (QGP). A new 
frontier in the study of QCD matter opened with the first 
collisions of 208 Pb ions in November 2010, at the Large 
Hadron Collider (LHC) at CERN. These collisions are 
expected to generate matter at unprecedented tempera- 
tures and energy densities in the laboratory. 
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The first step in characterizing the system produced 
in these collisions is the measurement of the charged- 
particle pseudo-rapidity density, which constrains the 
dominant particle production mechanisms and is essen- 
tial to estimate the initial energy density. The depen- 
dence of the charged-particle multiplicity density on en- 
ergy and system size reflects the interplay between hard 
parton-parton scattering processes and soft processes. 
Predictions of models that successfully describe parti- 
cle production at RHIC vary by a factor of two at the 
LHC [US- 

This Letter reports the measurement of the charged- 
particle pseudo-rapidity density produced in Pb-Pb col- 
lisions at the LHC, utilizing data taken with the ALICE 
detector [3] at a centre-of-mass energy per nucleon pair 
y / SNN = 2.76 TeV. The primary charged-particle density, 
dA^h/dry, in central (small impact parameter) Pb-Pb 
collisions is measured in the pseudo-rapidity interval 
\r]\ = | — lntan(#/2)| < 0.5, where 8 is the polar an- 
gle between the charged-particle direction and the beam 
axis (z). We define primary particles as prompt parti- 
cles produced in the collision, including decay products, 
except those from weak decays of strange particles. 

The present measurement extends the study of par- 
ticle densities in nucleus-nucleus collisions into the 
TeV regime. We make comparisons to model predic- 
tions [4rtl6j. and to previous measurements in nucleus- 
nucleus collisions at lower energies at the SPS and 
RHIC [T7H23] . as well as to pp and pp collisions over 
a wide energy range [2"6"H3"2"] . Our measurement provides 
new insight into particle production mechanisms in high 
energy nuclear collisions and enables more precise model 
predictions for a wide array of other observables in nu- 
clear collisions at the LHC. 

A detailed description of the ALICE experiment is 
given in Ref. [3]. Here, we briefly describe the detec- 
tor components used in this analysis. The Silicon Pixel 
Detector (SPD) is the innermost part of the Inner Track- 
ing System (ITS). It consists of two cylindrical layers of 
hybrid silicon pixel assemblies positioned at radial dis- 
tances of 3.9 and 7.6 cm from the beam line, with a total 
of 9.8 x 10 6 pixels of size 50 x 425 fim 2 , read out by 
1200 electronic chips. The SPD coverage for particles 
originating from the center of the detector is I77I < 2.0 
and \rj\ < 1.4 for the inner and outer layers, respec- 
tively. Each chip provides a fast signal if at least one 
of its pixels is hit. The signals from the 1200 chips are 
combined in a programmable logic unit which supplies a 
trigger signal. The fraction of SPD channels active dur- 
ing data taking was 70% for the inner and 78% for the 
outer layer. The VZERO detector consists of two arrays 
of 32 scintillator tiles placed at distances z = 3.3 m and 
z = —0.9 m from the nominal interaction point, cov- 
ering the full azimuth within 2.8 < r\ < 5.1 (VZERO- 
A) and -3.7 < 77 < -1.7 (VZERO-C). Both the ampli- 
tude and the time signal in each scintillator are recorded. 
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The VZERO time resolution is better than 1 ns, allow- 
ing discrimination of beam-beam collisions from back- 
ground events produced upstream of the experiment. 
The VZERO also provides a trigger signal. The Zero De- 
gree Calorimeters (ZDCs) measure the energy of specta- 
tors (non-interacting nucleons) in two identical detectors, 
located ±114 m from the interaction point. Each ZDC 
consists of two quartz fiber sampling calorimeters: a neu- 
tron calorimeter positioned between the two beam pipes 
downstream of the first machine dipole that separates 
the two charged particle beams, and a proton calorime- 
ter positioned externally to the outgoing beam pipe. The 
energy resolution at beam energy is estimated to be 11% 
for the neutron and 13% for the proton calorimeter, re- 
spectively. 

For the data analyzed, beams of four bunches, with 
about 10 7 Pb ions per bunch, collided at ^/snn = 2.76 
TeV, with an estimated luminosity of 5 x 10 23 cm _2 s -1 . 
The trigger was configured for high efficiency for hadronic 
events, requiring at least two out of the following three 
conditions: i) two pixel chips hit in the outer layer of the 
SPD, ii) a signal in VZERO- A, iii) a signal in VZERO- 
C. The threshold in the VZERO detector corresponds 
approximately to the energy deposition of a minimum 
ionizing particle. The luminous region had an r.m.s. 
width of 5.9 cm in the longitudinal direction and 50 /im 
in the transverse direction. The estimated luminosity 
corresponds to a hadronic collision rate of about 4 Hz. 
The observed rate was about 50 Hz, mainly due to elec- 
tromagnetically induced processes [33]. These processes 
have very large cross sections at LHC energies but gen- 
erate very low multiplicities and therefore do not con- 
tribute to the high particle multiplicities of interest for 
the present analysis. The trigger rate without beam was 
negligible and the rate in coincidence with bunches of 
only one beam was about 1 Hz. This beam background 
is eliminated from the triggered event sample using the 
VZERO timing information, as well as the correlation 
between the number of tracks reconstructed in the Time 
Projection Chamber (TPC) and the number of hits in 
the SPD. 

Offline event characterization utilizes global event ob- 
servables that are intrinsically correlated over different 
regions of phase space through the initial collision geom- 
etry. Figure [I] (upper) shows the measured correlation 
between the energy deposited in the ZDC and the sum 
of amplitudes in the VZERO detector. The VZERO re- 
sponse is proportional to the event multiplicity, and the 
ZDC energy to the number of non-interacting nucleons 
close to beam rapidity. As events become more central, 
with smaller impact parameter, they generate larger mul- 
tiplicity in VZERO and less energy forward in the ZDC. 
This behavior is understood based on collision geome- 
try and nuclear breakup [Mj. For small ZDC response 
the VZERO signal has two distinct values correspond- 
ing to peripheral and central collisions. However, the 
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FIG. 1. Upper panel: Correlation of ZDC and the VZERO 
response in hadronic collisions. Lower panel: Distribution of 
the sum of amplitudes in the VZERO scintillator tiles (black 
histogram); inset shows the low amplitude part of the distri- 
bution. The red line shows the fit of the Glauber calculation 
to the measurement. The fit was performed above the cut 
indicated in the inset, avoiding the region at low amplitudes 
dominated by the electromagnetic processes. The shaded area 
corresponds to the most central 5% of hadronic collisions. 



VZERO signal alone can be used to discriminate on cen- 
trality. Figure [T] (lower) shows the distribution of the 
VZERO amplitude for all triggered events after beam 
background removal. The distribution is fit using the 
Glauber model [35] to describe the collision geometry and 
a Negative Binomial Distribution (NBD) to describe par- 
ticle production [34] . In addition to the two parameters 
of the NBD, there is one free parameter that controls 
the power-law dependence of particle production on the 
number of participating nucleons (iV par t). To avoid the 
region contaminated by electromagnetic processes, which 
constitutes over 90% of the triggered events, the fit is 
restricted to the VZERO amplitude region above 150, 
where the trigger for hadronic collisions is fully efficient. 
The fraction of the hadronic cross section from the model 
fit corresponding to this cut, 87%, allows the determina- 
tion of the cross section percentile for any more-central 
VZERO cut by integrating the measured distribution. 
The most central 5% fraction of the hadronic cross sec- 
tion was determined in this way. 

The analysis is based on the VZERO event selection as 
described above. Among the triggered sample of about 
650000 events, 3615 events correspond to the most cen- 
tral 5% of the hadronic cross section, indicated by the 
shaded region in Fig. [I] (lower). The first step in the 
measurement of the charged particle multiplicity is the 
determination of the primary vertex position by correlat- 
ing hits in the two SPD layers. All events in the central 
sample are found to have a well-constrained primary ver- 



8 



tex. To minimize edge effects at the limit of the SPD 
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FIG. 2. Distribution of the azimuthal separation (Aip) of all 
candidate tracklets in data, simulation, and the background 
contribution that is estimated from the simulation. 

acceptance we have only used events with reconstructed 
vertex at |z v tx| < 7 cm, reducing the sample to 2711 
events. Tracklet candidates are formed using infor- 
mation on the position of the primary vertex and of hits 
on the SPD layers. A tracklet is defined by a pair of hits, 
one on each SPD layer. Using the reconstructed vertex as 
the origin, we calculate the differences in azimuthal (Aip, 
bending plane) and polar (AO, non-bending direction) 
angles for pairs of hits [3T]. Only hit combinations sat- 
isfying a selection on the sum of the squares of Aip and 
AO, each normalized to its estimated resolution (60 mrad 
in Aip and 25sin 2 mrad in AO), are selected as track- 
lets. If multiple tracklet candidates share a hit, only 
the combination with the smallest sum of squares of Aip 
and AO is kept. The cut imposed on Aip efficiently se- 
lects charged particles with transverse momentum (p t ) 
above 50 MeV/c. Particles below 50 MeV/c are mostly 
absorbed by material. 

The charged-particle pseudo- rapidity density d-/V c h/d?7 
is obtained from the number of tracklets within \rj\ < 0.5 
according to dN ch /dri = a X (1 - /3)dAtrackicts/d?7, where 
a is the correction factor for the acceptance and efficiency 
for a primary track to generate a tracklet and (3 is the 
probability to form a background tracklet from uncorre- 
cted hits. The corrections a and j3 are determined as a 
function of the z-position of the primary vertex and the 
pseudo-rapidity of the tracklet. The simulations used to 
calculate the corrections are based on the HIJING [3"6"] 
event generator and a GEANT3 [37. model of the detec- 
tor response. Three different methods have been used to 
estimate the combinatorial background. 

The main method to estimate the combinatorial back- 
ground j3 relies on the event simulation using a sample of 
events with similar multiplicities (SPD hits) as in the real 
data. In Fig. [2] the Ap distribution for candidate track- 



lets is compared for data and simulation. The distribu- 
tions are very similar, practically identical in the back- 
ground dominated tails. The second method is based 
on the injection of random background hits in the real 
event, in order to evaluate the probability of creating 
fake tracklets by combinatorics. In the third method 
events are modified by rotating hits in the inner SPD 
layer by 180° in ip, thereby destroying real correlations, 
but preserving global event features. In all cases, the ab- 
solute amount of combinatorial background is obtained 
by matching the tracklet and background distributions 
in the tails. For the main method, which ideally pro- 
vides both the shape and the normalization, an adjust- 
ment of 1% is needed to match the tails. The estimated 
combinatorial background is about 14%. In order to ac- 
count for the effect of the correlated background, the 
same background subtraction procedure is also applied 
to the simulation (i.e. without relying on the event gen- 
erator information). The correction for the acceptance 
and efficiency, a, is obtained by the ratio of the number 
of generated primary charged particles to the number of 
reconstructed tracklets after subtraction of the combina- 
torial background. In this way, a accounts for geomet- 
rical acceptance, detector and reconstruction efficiencies, 
contamination by weak decay products of strange parti- 
cles, conversions, secondary interactions and undetected 
particles below 50 MeV/c transverse momentum. The 
overall correction factor a varies slightly depending on 
vertex position and r/, and is about 2. 

We have considered the following sources of system- 
atic uncertainties: background subtraction estimated as 
2% by comparing the results of different methods; par- 
ticle composition estimated as 1% by changing the rela- 
tive abundances of protons, pions, kaons by a factor of 
two; contamination by weak decays estimated as 1% by 
changing the relative contribution of the yield of strange 
particles by a factor of two; \ow-p t extrapolation esti- 
mated as 2% by varying the amount of undetected par- 
ticles at low p t by a factor of two; event generator esti- 
mated as 2% by using HIJING [SS] with and without jet 
quenching, as well as DPMJET (38] for the corrections; 
centrality definition estimated as 3% by using an alter- 
native event selection based on the SPD hit multiplici- 
ties, and by varying the range of the Glauber model fit. 
All other sources of systematic errors considered (track- 
let cuts, vertex cuts, material budget, detector efficiency, 
background events) were found to be negligible. The 
total systematic errors amounts to 4.8%. Independent 
cross-checks performed using tracks reconstructed in the 
TPC and ITS yield results consistent within the system- 
atic uncertainty. 

In order to compare bulk particle production in dif- 
ferent collision systems and at different energies, and to 
compare with model calculations, the charged particle 
density is scaled by the number of participating nucle- 
ons, determined using the Glauber model fit described 
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FIG. 3. Charged particle pseudo-rapidity density per partic- 
ipant pair for central nucleus-nucleus [17H25] and non-single 
diffractive pp (pp) collisions 26-32), as a function of ^/snn- 
The solid lines oc s^n 5 an< A x Snn are superimposed on the 
heavy-ion and pp (pp) data, respectively. 



above (Fig. [TJ. The average number of participants for 
the 5% most central events is found to be (-/V part ) = 
381 with an r.m.s. of 18 and a systematic uncertainty of 
1%. The systematic uncertainty was obtained by vary- 
ing the parameters of the Glauber calculation within the 
experimental uncertainty and by ±8% around 64 mb for 
the nucleon-nucleon cross section, by using different fit 
ranges, and by comparing results obtained for different 
centrality variables (SPD hits, or combined use of the 
ZDC and VZERO signals). 

We measure a density of primary charged particles 
at mid-rapidity diV ch /dry = 1584 ± 4 (stat.) ± 76 
(sys.). Normalizing per participant pair, we obtain 
dAT ch /d7 ? /(0.5(iV paxt )) = 8.3 ± 0.4 (sys.) with negligi- 
ble statistical error. In Fig. |3j this value is compared 
to the measurements for Au-Au and Pb-Pb, and non- 
single diffractive (NSD) pp and pp collisions over a wide 
range of collision energies p"7lt32] . It is interesting to 
note that the energy dependence is steeper for heavy-ion 
collisions than for pp and pp collisions. For illustration, 
the curves oc s^sf and oc s^n 1 shown superimposed 
on the data. A significant increase, by a factor 2.2, in the 
pseudo-rapidity density is observed at ^/snn = 2.76 TeV 
for Pb-Pb compared to ^/snn = 0.2 TeV for Au-Au. 
The average multiplicity per participant pair for our cen- 
trality selection is found to be a factor 1.9 higher than 
that for pp and pp collisions at similar energies. 

Figure [4] compares the measured pseudo-rapidity den- 
sity to model calculations that describe RHIC measure- 
ments at y / SNN = 0.2 TeV, and for which predictions at 
^/snn = 2.76 TeV are available. Empirical extrapolation 
from lower energy data [1] significantly underpredicts the 
measurement. Perturbative QCD-inspired Monte Carlo 
event generators, based on the HIJING model tuned to 



7 TeV pp data without jet quenching [5], on the Dual 
Parton Model [6], or on the Ultrarelativistic Quantum 
Molecular Dynamics model [7] are consistent with the 
measurement. Models based on initial-state gluon den- 
sity saturation have a range of predictions depending on 
the specific implementation 8 -12 , and exhibit a varying 
level of agreement with the measurement. The prediction 
of a hybrid model based on hydrodynamics and satura- 
tion of final-state phase space of scattered partons [T3] 
is close to the measurement. A hydrodynamic model in 
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FIG. 4. Comparison of this measurement with model predic- 
tions. Dashed lines group similar theoretical approaches. 

which multiplicity is scaled from p+p collisions overpre- 
dicts the measurement |14j , while a model incorporating 
scaling based on Landau hydrodynamics underpredicts 
the measurement [15] . Finally, a calculation based on 
modified PYTHIA and hadronic rescattering [IB] under- 
predicts the measurement. 

In summary, we have measured the charged-particle 
pseudo-rapidity density at mid-rapidity in Pb-Pb colli- 
sions at y / SNN = 2.76 TeV, for the most central 5% frac- 
tion of the hadronic cross section. We find diV c h/d?7 = 
1584 ± 4 (stat.) ± 76 (sys.), corresponding to 8.3 ± 
0.4 (sys.) per participant pair. These values are signif- 
icantly larger than those measured at RHIC, and indi- 
cate a stronger energy dependence than measured in pp 
collisions. The result presented in this Letter provides 
an essential constraint for models describing high energy 
nucleus-nucleus collisions. 
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